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Abstract 


The paper treats the validity problem of the nonrelativistic Vlasov- 
Poisson equation in d > 2 dimensions. It is shown that the Vlasov- 
Poisson dynamics can be derived as a combined mean field and point- 
particle limit of an N-particle Coulomb system of extended charges. 
This requires a sufficiently fast convergence of the initial empirical 


distributions. If the electron radius decreases slower than N d ( d + 2 ). 
the corresponding initial configurations are typical. This result entails 
propagation of molecular chaos for the respective dynamics. 


1 Introduction 


We are interested in a microscopic derivation of the Vlasov-Poisson 
dynamics in d > 2 spatial dimensions. This is the system of equations 


dtf + p-V q f + (k*p t )-V p f = 0, 


(1) 


where k is the Coulomb kernel 




( 2 ) 


and 
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is the charge density induced by the distribution function f(t,p, q ) > 0, 
describing the density of particles with position q £ M. d and momentum 
p £ M d . Here, units are chosen such that all constants, in particular 
the mass and charge of the particles, are equal to 1. 

The Vlasov-Poisson equation provides an effective description of 
a collisionless plasma with electrostatic (a = +1) or gravitational 
(<r = —1) interactions. In the gravitational case, the equation is also 
known as Vlasov-Newton. 

Derivation of mean field equations. Kinetic equations of the 
Vlasov type are usually conceived of as mean field equations , effec¬ 
tive descriptions of many-particle systems, in which the IV-particle 
interactions are approximated by an “average” effect, determining an 
autonomous time-evolution for the distribution function /. 

Classical results, dealing with simplified models with Lipschitz- 
continuous forces, prove a statement of the following kind: If for 

an initial microscopic configuration Z = (qi,Pi)i=i...,N the empiri- 

N 

cal distribution Pq [Z] = jf d qi 5 Pi approximates a continuous den- 

i —1 

sity /o, then, at time t > 0, the time-evolved distribution p^ = 
N 

~k bqi(t)^Pi{t) approximates ft, where ft is a solution of the corre- 
i=l 

sponding Vlasov equation of the form ©• Formally, the approximation 
is understood in terms of weak convergence of probability measures, 
quantified by an appropriate metric. (Neunzert and Wick, 1974 [19], 
Braun and Hepp, 1977 @], Dobrushin, 1979 |6|; see also fTSl 124] .) 

Results for singular forces. For singular forces - up to but not 
including the Coulomb case - similar results could be proven only 
recently. Hauray and Jabin (2013) [10] treat force kernels bounded 
as \k(q)\ < with a < d — 1. For 1 < a < d — 1 they require 

an IV-dependent cut-off which can be chosen as small as N~ l / 2d for 
a yZ d — 1, while for a < 1 they are able to perform the mean field 
limit without cut-off. Pickl and Boers (2015) improve the cut-off for 
singularities near the Coulomb case to N~ l / d [lj. 

The microscopic justification of the Vlasov-Poisson equation, cor¬ 
responding to the case a = d— 1, has been an open problem, so far. In 
this paper, we propose a particle approximation by extended charges 
with IV-dependent radius that can decrease as fast as N d( - d +v +c . The 
proof is based on a stability result of Loeper H9 and an anisotropic 
variant of the Wasserstein distance. An alternative proof, general¬ 
ization the methods introduced in |T], is simultaneously presented in 
HU. 
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While the proof in m requires a much smaller cut-off of order 
N~d +e , the method presented here allows for better rates of conver¬ 
gence. Moreover, the microscopic regularization proposed here can be 
understood as a nonrelativistic analogue of the rigid charges model 
that was used by Golse (2012) to perform the mean field limit for a 
regularized version of the Vlasov-Maxwell dynamics [8]. Our discus¬ 
sion might thus also be interesting in view of a possible generalization 
to the relativistic Vlasov-Maxwell system. 

Convergence in law and molecular chaos. What all recent results 
with singular forces have in common, is that they are probabilistic in 
the sense that the mean held limit can be performed for typical initial 
conditions. In other words, the microscopic density ^ converges in 
law to the constant variable ft, which is given as the solution of the 
corresponding mean held equation. By a well-known result in prob¬ 
ability theory (e.g. j 13j , [9], [25] Prop.2.2 ], |T7]). this is equivalent 
to molecular chaos in the following sense: If at time t = 0 the parti¬ 
cles are identically and independently distributed with law /o and if 
the corresponding product measure Fq = <8> W /o on M 6JV evolves with 
the microscopic iV-particle how, then, for times t > 0, it holds that 
F t N = ~ <g) N ft, where the approximation is understood in 

terms of convergence of marginals. That is, writing Xi = ( qi,Pi ), we 
consider the fc-particle marginal 

( k) F t N (xx k ) := j F t N (Z) dx k+1 ...dx N . 

Then (^F^ converges weakly to <8> fc /o as N —>• oo for all k £ N. 

Note that the probabilistic nature of these statements is in contrast to 
earlier results in [3] and |6j, which are, in effect, deterministic, allow¬ 
ing arbitrary sequences of initial configurations. One reason is that for 
singular forces, there exist “bad” initial conditions leading to cluster¬ 
ing of particles and significant deviations from the typical mean held 
behavior. 

Comparison of recent result, open problems The strategy em¬ 
ployed in [10], as well as in the present paper, is thus to impose addi¬ 
tional constraints on the initial conditions, subsequently showing that 
these constraints are satished with probability 1 in the limit IV —>• oo. 
In [10], the respective bounds are imposed on the concentration of 
particles at t = 0, while in our proof, the probabilistic element enters 
through the requirement of a sufficiently fast convergence of the initial 
microscopic distribution. In any case, these assumptions assure that 
the initial configuration is “well-placed”, so to speak, preventing, in 


3 


particular, a blow-up of the microscopic dynamics. 

One of the key innovations in our proof is that a regularization 
is applied on the level of the charge density, which allows us to work 
with bounded densities rather than Dirac masses. The L°°-norm of the 
microscopic charge density is controlled by propagating a respective 
Wasserstein bound (Lemma 14.31) . Similar estimates can be used to 
carry over stronger regularity properties from the Vlasov density to the 
microscopic density. While this method is rather simple, it requires a 
relatively large cut-off of order ~ N d(d+2) +c for some e > 0. Moreover, 
there is no immediate connection between the width of the cut-off and 
the strength of the singularity. This is in contrast to the situation 
in m or [T], where the lower bound on the required regularization 
decreases with a. 

Most notably, Hauray and Jabin, 2013 [10] . are able to prove molec¬ 
ular chaos for weak singularities (a < 1) with no cut-off at all, while 
the method proposed here requires in any case a regularization (smear¬ 
ing of the charges). However, we emphasize that the results in [1U] do 
not include the Coulomb case a = d — 1, which is the main focus of 
our paper. Furthermore, our result applies also in dimension 2, while 
the assumptions required in pLO] are no longer generic in that case. 

The method introduced in fl[ and extended in [13] is designed for 
stochastic initial conditions, thus aiming directly at a typicality result. 
Rather than controlling the difference between f N and fi N in some 
weak metric, one considers a stochastic process of the form E(| 

$co( Z )l°o), where are the N-particle flows generated by the 

mean field dynamics and the microscopic dynamics, respectively. The 
corresponding proof allows the cut-off to decrease as fast as N~ l / d+e , 
i.e. (almost) as fast as the typical distance between a particle and its 
nearest neighbour. 

Whether mean field results for strongly singular forces - approach¬ 
ing or even including the Coulomb case - can be obtained with no 
cut-off at all is an open question. Concerning related problems, we be¬ 
lieve that some of the methods presented here could be generalized to 
the relativistic Vlasov-Maxwell dynamics and leave this to be treated 
in a future paper. 


2 The microscopic model 

As the force kernel considered here is strongly singular at the origin, 
we will require a regularization on the microscopic level. We shall 
consider as a microscopic model the dynamics of smeared (extended) 
charges with Coulomb interactions. The cut-off parameter rjsr thus has 
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a straight-forward physical interpretation as a finite electron radius. 
In the relativistic case, an analogous model of rigid charges (without 
collisions and rotations) was used by Golse to derive a regularized 
version of the Vlasov-Maxwell system [5] (c.f. also Rein, 2004 [ST]). 

While the smearing of charges is a natural way to regularize point- 
interactions, the cut-off thus introduced is a technical necessity rather 
than a realistic physical model. In particular, the ^dependence of 
the electron radius might seem strange from a physical point of view, 
though similar regularizations are commonly used in numerical simula¬ 
tions. Also note that in our proof, the radius has to be chosen so large 
that a great number of particles will typically overlap. Intuitively, the 
combined limit N —>• oo, rjv —>• 0 describes a regime where a large 
number of smeared electrons blurs into a continuous charge cloud. 


As before, let 

k:R d ->R d , q^cr-^-, 

\q\ d 

denote the Coulomb kernel. That is, if T : W l —>• M is a solution of 
Poisson’s equation 

AT = cp, lim T(g) = 0 

|<?|—H-oo 

in the sense of 

^(?) = [ i- -7TcT2 P(^)dV, d > 3, 

J \q- q'\ a z 

or 

T(g) = —a j ln(g — q')p(q') d q', for d = 2 

then 

-VT(qr) = k * p{q) =cr J J _ V p(q') d d q'. 

We consider a system of N charges, smeared out by a smooth, non¬ 
negative, spherically symmetric form-factor x € C , Q°(M rf ). For simplic¬ 
ity, we shall assume that x satisfies: 

i) supp(x) A B(1;0) = {i 6 M d : ||a:|| < 1}. 
h) llxlloo = SUPa.gR x(x) = 1. 
iii) llxlli = f x{ x )dx = 1 . 

We call a sequence (r/v)iVe n of positive real numbers a rescaling se¬ 
quence if it is monotonously decreasing with r\ = 1 and lim = 0. 

TV—>-oo 

Given such a rescaling sequence, we define a rescaled form-factor as 

X N ( X ) : = \x(—\ N € N. (4) 

r d N \r N J 
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The configuration of the microscopic system is given by Z(t) = (<7j(t),Pi(t))i<i<jv, 
where qfit) is the center of mass of particle i, and pfit) the correspond¬ 
ing momentum at time t. The equations of motion in the so called 
mean field scaling read: 


with 


4i(t) =Pi{t ) 

Pi{t) = K N (q. i ;q 1 ,...,q N ) 


(5) 


K N (q%,qu-,qN) ■= j J x N (qj-y)k(z-y)x N (qi-z)d d yd d z. 

3 =1 

( 6 ) 


The IV-particle force © can be rewritten in the following way: Given 

N 

the microscopic density p,™ = jj 4 qi (t)4 Pi (t), one checks that 

i =1 


K N (-,q i, Qn) = X N * k * X N * p[Pt\ =■ k * p[$]. 


where we introduce the notation 

P ■= X N *x P (7) 

for (p, a measure or measurable function on W l and * denoting the 
convolution with respect to the space variable. 

Except for the scaling-factor IV -1 , these equations describe the 
regular Coulomb dynamics for smeared charges with form-factor x N ■ 
They can be understood as the nonrelativistic limit of a Maxwell- 
Lorentz system of rigid charges (also known as the Abraham model, 
c.f. [23] Chs. 2,13]). The double-convolution results from the fact that 
the charge enters the interaction-term quadratically; In other words, 
the charges acting and the charge being acted upon are both smeared 
out. Note that this system is Hamiltonian for 

jv i 

H(q t ,Pi) = Y, 2 P * + 2n12 / / *(P ~ qi \~ °\d-i X( z ~ Qj) d y 
i= 1 ^ i,j J J \ y\ 

and thus conserves total energy. Note also that this Hamiltonian in¬ 
cludes self-interactions. 
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2.1 The regularized Vlasov-Poisson equation 

For the microscopic model described above, we introduce a correspond¬ 
ing mean field equation: 


dtf+p-V q f + k N [pt]-V p f = 0 , 

k N [pt](q) ■■= x N *k* X N * pt(q), ( 8 ) 

Pt(q) = p[ft\{q) = j f(t,q,p)d d p 

we call this the regularized Vlasov-Poisson system with cut-off parame¬ 
ter r/v- For N —>• oo, the form-factor \ N approximates a delta-measure 
in the sense of distributions and © formally reduces to the Vlasov- 
Poisson system ©• 


2.2 The method of characteristics 

Let v = (vt)te[o,T) a continuous family of probability measures on M. d X 
for T e M + U{+oo}. Let pt\y\(q) = f v(q,p) d d p the induced charge 
distribution on W 1 . We denote by gf t s = (Q u (t,s,qo,po),P l '(t,s,qo,po)) 
the one-particle flow on M. d X solving: 

'itQ = p 

{ it p = X N * k * X N * p(Q) 

Q{s,s,q 0 ,Po) = qo 

P(s,s,q 0 ,p 0 ) =Po 

This flow exists and is well-defined since the vector-field is Lipschitz 
for all N. If f N (t,q,p ) is a solution of ©, it is now straight-forward 
to check that 

ft N = vC#f f, VO < s < t < T. (10) 

Here, <^(-)#/ denotes the image-measure of / under (p, defined by 
ip#f(A) = for any Borel set A C M 6 . 

Conversely, if ft is a fixed-point of (vf) —>• it is a solution 

of (J8l) with initial datum fo. In particular, one observes that Z(t) = 

N 

(qi(t),Pi{t))i=l,..,N is a solution of © if and only if p N [Z(t)] = jj S qi {t)$Pi(t) 

2—1 

solves © in the sense of distributions. Basically, our aim is thus to 
show that this relation carries over to the limit N —>• oo. 

For the (unregularized) Vlasov-Poisson equation, the corresponding 
vector-field is not Lipschitz, in general. However, if we assume the 


7 



existence of a solution ft with p £ L°°([0, T] x M. d ), the mean held force 
k* pt does satisfy a Log-Lip bound of the form \k * pt(x) — k* pt(y)\ < 
C\x — y||log(|x — 2 /|)| (for \x — y\ < b, let’s say). This is sufficient to 
ensure the existence of a characteristic how = ( Qt,s,Pt,s ) solving 

dtQt,s = Pt,s 
< d t.Pt,s = k * p[ft\(Qt,s) 

Q(s,s,q 0 l Po) = qo 
P(s,s,q 0 ,p 0 ) =p 0 


such that ft = ipt,sibfsi f° r all 0 < s < t < T. 

2.3 Existence of Solutions 

For the regularized mean held equations (j5J), all forces are Lipschitz 
continuous and the solution theory is fairly standard, see e.g. Earn. 
For the actual Vlasov-Poisson system, the issue is more subtle. For¬ 
tunately, in the physically most relevant, 3-dimensional case, we can 
rely on various results, establishing global existence and uniqueness of 
(classical) solutions under reasonable conditions on the initial /q (Pfaf- 
felmoser, 1990 |20| . Schaeffer, 1991 [22], Lions and Perthame, 1991 [15]. 
Horst, 1993 [12]). The situation is similar in the 2-dimensional case, 
treated in Ukai and Okabe, 1978, [25] an( i Wollmann, 1980 [28] . 

For the rest of the paper, we shall work under the following assumption: 

Assumption 2.1. For /o £ L 1 nL°°(W i xW 1 -, Rg") there exists a T* > 0 
such that the Vlasov-Poisson system (1-3) has a unique solution ft on 
[0, T*) with /(0, •,•) = /o- Moreover, as we consider the sequence of 
solutions to the regularized equations (]8|, the charge density remains 
bounded uniformly in N and t, i.e. 3C*o < +oo such that 

||pL/f]||oo<0), vt <T* VA£NU{+oo}, (12) 

where, with a slight abuse of notation, ff° := ft. 

In fact, given a bounded charge density, uniqueness of the solution 
(in the set of bounded, positive measures) is proven in Loeper, 2006 
[16] . Moreover, it is well known that as long as the charge density is 
bounded, solutions with smooth initial data remain smooth (see e.g. 
in [Ft]). 

In dimension d = 3, the existence result of Lions and Perthame [T5] 
ensures that the above assumption is satisfied for a relatively large class 
of initial data and T* = +oo. 
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Theorem 2.2 (Lions and Perthame). 

Let /o >0,/o 6 L 1 (M 3 x R 3 ) n L°°(R 3 x R 3 ) satisfy 

J\p\ m fo(q,p)dqdp < Too (13) 

for all m < mo and some mo > 3. 

1) Then, the Vlasov-Poisson system defined by equations (1-3) has 
a continuous, bounded solution f(t , ■,•) € C(R + ;L r (R 3 X R 3 ) D 
L°°(R + ; L°° (R 3 x R 3 )), 1 < r < oo, satisfying 

sup / \p\ m f(t,q,p)dpdp < Too, (14) 

te[o,T] J 

for all T < oo, m < mo- 

2) If, in fact, mo > 6 and we assume that fo satisfies 

supess{/ 0 (g / +pt,p') : \q - q'\ < Rt 2 , \p-p'\ < Rt} 

€ L°°((0,T) x Rg;L 1 (R 3 )) 

for all R > 0 and T > 0, then there exists C > 0 such that 

l|pL/f]||oo<C', Vt>0VJV€NU{oo}. (16) 

Note that Lions and Perthame state m only for f t , though they re¬ 
mark (and it is straightforward to check) that the proof actually yields 
an upper bound on the charge densities p[/( v ] as one considers a se¬ 
quence of regularized time-evolutions as, for instance, in (J8]). 

In higher dimensions, where blow-up might occur, there exists at least 
some T* > 0, depending only on fo, such that (11511 is satisfied, if one 
assumes that fo has compact support. This is ensured by the following 
lemma. 

Lemma 2.3 (Local existence of solutions). Let fo € L 1 nL°°(R d x R 3 ) 
with compact support and f a (local) solution of Jl]) with f\t=o= fo- 
Let 

D(t) := sup{|g| :3p€R d : f(t,q,p) 0} (17) 

R(t) := sup { \p\ : 3q € R d : f(t,q,p) 0} (18) 
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the diameter of the support in the q-, respectively p-coordinates. Then 
there exists a constant C > 0 such that 

t 

D(t) < D(0) + J R(s)ds (19) 

o 

t 

R(t)<R(0) + C\\f 0 \\ oo \\f 0 \\ 1 1 /d J R d ~\s)ds. (20) 

o 

These estimates hold independent of N as we consider the sequence 
f N of solutions to the regularized equation ((HJ) with f N \t=o= /o- 

3 Statement of the results 

Our approximation result for the Vlasov-Poisson dynamics is formu¬ 
lated in terms of (modified) Wasserstein distances. In the context 
of kinetic equations, the Wasserstein distance was introduced by Do- 
brushin [6j. Here, we shall briefly recall the definition and some basic 
properties. For further details, we refer the reader to the book of 
Villani |27J Ch. 6]. 

Definition 3.1. Let V(M. k ) the set of probability measures on WL k 
(equipped with its Borel algebra). For given /i, v E V(M. k ) let n(/x, v) 
be the set of all probability measures x with marginal /r and v 
respectively. 

For p E [1, oo) we define the Wasserstein distance of order p by 

W p {ii,v):= inf ( [ \x - y\ p dvr(x, y) ) /P . (21) 

R fc xR fc 

Convergence in Wasserstein distance implies, in particular, weak con¬ 
vergence in V(M. k ), i.e. 

J $(x) dpL n {x) —> J <L(a;) d/i(.x), n —>• oo, 

for all bonded, continuous functions <F. Moreover, convergence in W p 
implies convergence of the first p moments. W p satisfies all properties 
of a metric on V(M. k ), except that it may take the value +oo. 

An important result is the Kantorovich-Rubinstein duality: 

W p(TW) = sup{ J $i(x)dp(x) - j <5> 2 {y) dv(y) : 

($ 1 , $ 2 ) € L l {n) x L l (y), $i(y) - § 2 {x) <\x- y| p |. 

( 22 ) 
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A particularly useful case is the first Wassersteiu distance, for which 
the problem reduces further to 


W\{pL,v) = sup | / $(x)d/x(x)— / $(x) dv(x) 


where ||$||ij p := sup , to be compared with the bounded Lip- 

x^y x y 

schitz distance 

dBL{p, v) = sup| J $(x) dfj,(x) - J $>(x) du(x ); ||$||z,i P , ||$||oo < l}- 

We can now state our precise results in the following theorems. 

Proposition 3.2 (Deterministic Result). Let /o £ L 1 n L°°(R rf X 
R 1d ), / > 0. Let (rjv )NeN be a rescaling sequence and the unique so¬ 
lution of the regularized Vlasov-Poisson equation ([8| with f N { 0, -, -) = 
/o- Assume that on [0,T] the sequence {/n)n satisfies the uniform 
bound m on the induced charge densities. Suppose we have a se¬ 
quence of initial conditions Z £ R 6Ar such that 

Jim r- {1+ l +e) W 2 (^[Z]J 0 ) = 0 (23) 

for some e > 0. Then it holds that 

lim r“ {1+|) W 2 (/rf [Z\, f t N ) = 0, VO < t < T. (24) 


Since we will also show that W-fifl*, ft) = o(r^f e ) (Proposition 14. 6|) 
this establishes a particle approximation of the Vlasov-Poisson equa¬ 
tion for initial conditions satisfying fl23l) . 

Theorem 3.3 (Typicality Result). Let f o € L°°( R^xR 6 *) a probability 
measure such that the Vlasov-Poisson equation (HD has a uniqiLe solu¬ 
tion on [0, T*),T* £ R + U {+oo} with /(0, ■, •) = /o- Assume that the 
sequence (/iv)jv of solutions to the regularized Vlasov-Poisson equation 
Q with the same initial data satisfies the uniform bound m on the 
induced charge densities. Assume, in addition, that there exists k > 4 
such that 

M k (fo) := J (|<?| + \p\) k fo(q,p)<iqdp < +oo. (25) 

Suppose that r^ > N~ s with 

6 = d{2 + d + 2efi 6 > °' 
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Then there exist constants C\, C 2 , C 3 such that for all T < T* and N 
large enough that rjy < exp[—( 2 ( ' 1 J 1+1 ) 2 ] it holds that 


sup W 2 (p?[Z\J t ) > r]fi e < C 2 (e~ C3Ne + IV 1 "^), (26) 

4e[o,T] J 


where the probability Po defined in terms of the product measure 
fo on (M rf X M. d ) N . The constant C\ depends on d,x and Cq as in 
m , while C 2 , C 3 depend on d, k and 


Remark 3.4. 

1. In dimension 3, the necessary cut-off is of order N~ s with S < 

2. If the finite moment condition (1251) is replaced by the assumption 
of a finite exponential moment J e 7 ^”d/o(x), the rate of conver¬ 
gence becomes exponential, as well. This holds, in particular, for 
compactly supported /q . 


3.1 Sketch of the Proof 

We give here a brief sketch of our derivation and the central concepts 
and ideas on which it is based. 

1. To control the distance between microscopic density and mean 
field density, we introduce a variant WJfi of the second Wasser- 
stein distance defined with respect to the iV-dependent metric: 

d N {(qi,Pi), ( d 2 ,P 2 )) := (1 V \/|log(r/v)|) |gi - q 2 \ + \pi~P 2 \- 
where a V b := max{a, b}. 

2. We use an estimate from Loeper’s proof of uniqueness of weak 
solutions with bounded density [l 6 j to control the L 2 -norrn of 
the difference between mean field force and microsocpic force in 
terms of the quadratic Wasserstein distance. 

3. The regularization yields a Lipschitz bound on the microscopic 
force that diverges logarithmically with N. In terms of the mod¬ 
ified Wasserstein distance, this leads to a Gronwall estimate of 
the form 

< Cy/\\og(r N )\W 2 N (p?,f t N )- 

4. The previous bounds can be applied if the (smeared) microscopic 
charge density p^ = x N * p[pt\ remains bounded uniformly in N. 
We show that this can be assured as long as W 2 (pf [Z\, f^) = 
0{r~ {1+d,2) ). Given a sufficiently fast rate of convergence at 
t = 0, i.e. assumption (1231) . we conclude with 3. that this bound 
propagates. 
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5. It remains to check that the constraints so imposed on the ini¬ 
tial data are satisfied for typical Z % if the initial configuration 
is chosen randomly according to the product law <S> N fo- This is 
achieved with a recent large deviation estimate found by Fournier 
and Guilin [Tj. This estimate also sets the upper bound on the 
rate at which r/v can go to zero in the limit N —>• oo. 

4 A Gronwall-type Argument 

Our mean held limit is based on the following stability result by Loeper 
[161 Thm. 2.9], which is proven by methods from the theory of optimal 
transportation. 

Proposition 4.1 (Loeper). Let k the Coulomb kernel and pi, P 2 € 
L 1 (M. d ) D L°° (M. d ) two (probability) densities. Then 

\\k * pi - k * P 2 W 2 < [max{||pi|| 0O , ||p 2 ||oo}] 1/2 ^2(pi,/32)- (27) 

Moreover, we require the following estimates on the mean held force: 

Lemma 4.2. Let k as before and p € L 1 (M d ) 0 L°°(R d ). Then it holds 
that 

1 ) Ilfc^lU^I^IMoo + Mr. 

ii) ||*" * k* p\\ Lip < C L { 1 V |log(rjv)|) (||p||i + IIpIIoo) 

where we use again the notation a V 6 := max{a, b}. |S' c( “ 1 | denotes 
the area of the unit sphere and Cl is a constant depending on *. 


Proof, i) For the hrst inequality, we compute 


II A: * p||oo < 

J k(y)p t (x - y)d d y 

+ 

(X) 

J k(y)p( x - y) d d y 


|2/|<! 


\y\> 1 


< HHIco J j^p= r d' i » + w I = |s d - 1 |W«, + Wi 

?/l<1 


ii)We split the expression as 

II v(x *k* p)|l < II V(x * ^u 4+1 * p)L +1| v(x * K<rr *p>\L 

< ll.^lli II * p|L + IIVx-IL IlfcU^illr ||p|| 
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Now, we have: 


Vk \x>r d+1 *p( x ) 
—' N 


< / 1 -tffc-vWv 

\y\>ri +1 

< j ^A*-y)i d v+ / j^(* 

r-^ H1 <|j/|< 1 M >:L 

< (d + 1)C IIpIIoo log(r^) + ||p||i. 


y)d d y 


Furthermore: 

HVx^lloo = r- (d+1) ||Vx||oo 

and 

ll*U«IL= / j^d^ = i sArT- 

\y\<r d N +1 

Putting everything together, the statement follows. 


□ 


For the continuous solutions to the (regularized) Vlasov-Poisson 
equation, the corresponding charge-densities pt = p[ft*] are bounded 
by assumption. The challenge is to provide a bound on the microscopic 
charge density that holds uniformly in N. i.e. as the electron radius 
decreases and the forces become more singular. The idea is to show 
that as long as and are close as probability measures, the L°°- 
norm of p[f^] provides a bound on the L°°-norm of p[p^]- A simple 
such estimate was obtained in [3j Prop. 2.1] for the first Wasserstein 
distance. In view of the general Kantorovich-Rubinstein duality, we 
generalize this result to Wasserstein distances of higher order. 


Lemma 4.3. Let pi,p 2 two probability measures on W l and p 2 € 
L°°(R d ). Then: 


l|pi||oo<|B d (2)|||p 2 || co + r- (p+,i) W;(p 1 ,p 2 ), (28) 

where B fi (2) C M. d is the d-dimensional ball with radius 2. 

Proof. For any integrable function 4*, we consider the c-conjugate 


<J> c (y) := sup{<P(x) — \x — y\ p } 

X 

This is the smallest function satisfying 4> c (y) > <F(y) and <h(x) — 
& c (y) < \x — y \ p , Vx, y € M d . 
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Now, we write 


Pi(x) 


~(d+p) 
' N 


r d N P X N (x ~ y)pi(y)dy - j(r d + p X N (x - ■)) c (z)p 1 (z) dz 
+ f ( r< N P X N (x~-)) c (z)pi(z)dz 


By the Kantorovich duality theorem (1221) we have 
J r d jf P x N {x-y) pi{y)dy-J (r d + p X N (x— )) c (z) p 2 {z)dz < W p (pi, p 2 ). 
It remains to estimate 


J {r d N P X N {x - -)) c (2) P2(z) dz. 


Recalling that Hx^Hoo = x N d , we find 

(rfr V x N (x--)) c (z) = sup {r d ^ p X N (x-y)-\y-z\ p } < rJ^Hx^lloo = r p N . 

y& R 3 

Moreover, we observe that 

supp (r < ^~ p X N (x — -)) c C B(2r7v; x) : = {z € M 3 : \z — x\ < 2r at}, (29) 

since \z — x\ > 2implies X N (x — y) = 0, unless \y — z\ > rjv- But 
then: r^~ p X N (x — y) — \y — z\ p < r d j^ p rff — r p N = 0. Hence, 

J{r d N P X N {x ~ -)) c {z)p 2 {z)dz < 11/?211oo |B(2r7v; x)\ < 2 d |B d (l)| ||p 2 ||oo r d + p 
In total, we find 

llpilloo < r~ {p+d) W p ( Pl ,p 2 ) + \B d {2)\\\p 2 \\oo 


as announced. □ 

We shall apply the previous Lemma to p\ := p[p^(Z)\ and p 2 := p[f™] 
using \\p[f t N ]\\ < C p and W 2 {p[p?(Z)\, p[f t N }) < W 2 (p?(Z)J t N ) to get 
a bound on the (smeared) microscopic charge density. 

Finally, we need the following inequalities for the smeared densities. 

Lemma 4.4. 

Let x € C^°(M rf ), (?’Ar)iv a rescaling sequence and X N the rescaled 
form-factor as defined in ©• Let p, v E V(M. d ) and u = X N * x v etc. 
Then we have for 1 < p < oo: 

i) W p {i>,v) < r N 
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a) w p (ji, v) < w p (fx,u). 

Proof, i) Define ir(x, y) := v(x)x N (x—y) and observe that f dx 7 r (x,y) = 
P(y), f dy7r(x,y ) = v{x), hence it € II(i>, v). it has support in 
{\x — y\ < rjy}. Thus, we conclude 

W p {u, v)= inf f [ \x — y\ p dir'{x, y) \ 

R d xR d 

\ 1 /P 

x - y\ p dir{x,y) J <r N . 

R d xR d 



ii) In view of the Kantorovich duality (1221) . we find for ($ 1 ,^ 2 ) € 

L 1 (/r) x T 1 (^) with <J>i(y) — & 2 (x) < \x — y\ p : 

J $i(x) dp,(x) - J $ 2 ( 2 /) dz%) = J (x * $i)(a;) dp,(x) - J (x * $ 2 )(y) <My) 

But * $1 and \ * $2 also satisfy 

\x * $i(x) -X* ® 2 (y)\ = J x(z)$i(x - z) dz - J x(z)$ 2 (y - Z) dz 
< J x( z ) |$i(aJ - z) - $ 2 ( 2 / - z)\dz < J x(z) \x - y\ p dz = \x - y\ p . 

Hence, we have 

j d*! dft — j ^>2 d v = j $1 dfj, — J 4*2 < W p (n, u) 

and taking the supremum over all (<hi, < 1 * 2 ) yields the desired inequality. 

□ 


4.1 Modified Wasserstein distance 

As we want to establish a Gronwall inequality for the distance between 
empirical density and Vlasov density, we aim for a bound of the form: 

distfii^At, /t+ At ) - dist(/if, f t N ) oc dist(/rf, if) At + o(At). 

The choice of a metric, giving precise meaning to dist(/rf, /f), is thus 
a balancing act. While a stronger metric is, in general, more difficult 
to control, it also yields stronger bounds as it appears on the right 
hand side of the Gronwall estimate. 

If we compare the characteristic flow of the mean held dynamics 
with the how corresponding to the “true”, i.e. microscopic, dynamics, 
the growth in the spatial distance is trivially bounded by the distance 
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of the respective momenta. The only problem lies in controlling fluctu¬ 
ations in the force, i.e. the growth of the distance in momentum space. 
The idea, first employed in |[L4], is thus to be more rigid on deviations 
in the (/-coordinates, weighing them with an appropriate ^-dependent 
factor, and use this to obtain better control on the forces. 

Definition 4.5. Let ( r/v)iVeN be a rescaling sequence. On x W l 
we introduce the (^-dependent) metric: 

d N {(qi,Pi), ( Q2,P2 )) := (1 V v /|k ) g(rjv)|) \ qi - q 2 \ + \ Pl - p 2 \. (30) 

Now let Wp{-,-) be the p’th Wasserstein metric with respect to d N , 
i.e.: 

Wp(/J,,v):= inf ( f d N (x,y) p dn(x,y)] /P . (31) 

y rren^u) V J J 

R d xR d 

Note that W p (n, v) < W™(y, v)<(lV A/|log(rjv)|) W p (y, v), V/z, v € 
V(R d x M d ). Finally, we define 

W*{p,u) :=min|l, r JV (1+z) W?(y, z')}- (32) 

Obviously, convergence with respect to W* is much stronger than con¬ 
vergence with respect to W 2 . Concretely, we have for any sequence 
(^n)ngn and v € V(M d x M rf ): 

W*(v N ,v) ->■ 0 => WtiyN, v) = o(r^ +2 ). 


4.2 Deterministic result 


We now come to the central part of our argument: 

Proof of Proposition l3T2l Let IV € N and ttq € Let 

ipt = (Qt ■ P( l ) and ip{ = ( Q{,P ’/) the flow induced by the char¬ 
acteristic equation (j9j) for yf and fl' T . respectively. For any t € 
[0,T], T < T*, define the (^-dependent) measure 7 q on M. 6N X M. GN 
by 7 r t = Then n t € II (y?,f t ), Vt G [0,T], We set 


m ■■= 


/ 


d N (x,y) 2 d? T t (x,y) 



J ((IV Vllog^Tv)!) I® 1 

R 6 xR 6 


2 / 1 ! + \x 2 - y 2 


0 


2 


d7 T t (x,y) 



J ((1 v Vliog^Tv)!) | Qt(x) - Q{(y)\ + \ p t (x) - P/{y)\) 

R 6 xR 6 


2 


d7T 0 (x,y) 


1/2 
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Note that , f t N ) < D(t ) for any ttq € II(/o,/o). Now we con¬ 

sider: 

D*(t) :=mm{l, r~ {1+ ^ D(t)}. (33) 

Obviously, ^ D*(t) < 0 whenever D(t) > r^~ 2 since D*(t ) is already 
maximal. For D(t ) < 2 , we compute: 

- 

2 J ((1 V Vltogfav)!) |Qt*(x) - Q{(y) | + |ff (®) - P/(y)|)- 

((IV Vliog^Jv)!) - P/(y)| + |fc * Pt(Qt ( x )) “ * * p/(< 3/(2/))|) dvr 0 (x, y). 

The interesting term to control is the interaction term 

\k*Pt {•Qti x )) - £ * pf (Qt (y))| 

< - k*ft(Q{(y))\ (34) 

+ \k*Pt (Qtiv )) -k*pl ( Qt(y))\ ( 35 ) 

We begin with (1341) and hnd with Lemma 14.21 

\ k * Pt ( Qt( x )) ~~ k * Pt (Qtiv)) I , 3g * 

< c L (i v |iog(rjv)|)(i + HpKHoo) ~ <3f(y)| 

Hence, we have 

^ 2 (t)< Ji(t) + J 2 (t) (37) 

with 


Ji(t) := 2 J d7r 0 (x, y) ((lVy / jlog(r/v)j) \Q^(x)—Q{(y)\+\Pf l {x)~p/ (y)|) • 

((lV v / |log(r JV )|)|P/"(x)-P/ (y)|+C' L (lV|log(r iV )|)(l+||p( 1 ||oo)|Qi 1 (a:)-(5/(y)|) 

(38) 


Mt) := 21((1 V Vliog(nv)l) IQ^(*) - Qf(y)l + |if (*) - ^/(y)l)- 

* fit (Qt (y))-k* p{ (Q{(y)) I dvr 0 (x, y) 

(39) 

Now we observe that 

■M*) < <^(1 V |log(r’jv)|)(1 + |K||oo)£> 2 (t), (40) 
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while for the second term, we find with Holders inequality 

■h(t) < 

2 /((IV V|log(r-jv)|) | Qt(x) ~ Q{(y) I + lit*(*) - H 

\k*pt (Q{(v)) - ~ k * Pt {Qt\y ))\ 2 
We identify (HIT) as D(t), while for (TLZ1) we get 


d7r 0 (.T,y) 

(41) 

(42) 


k* (Pt ~ Pt){Qt(y))\ Z d? T 0 (x,y) 


1/2 


k * (.Pt ~ Pt) ( Qo(y ))r d7r t(x,y) 


1/2 


< 


k*ff -k* p{j 2 fit , y) d 2d y) 
k*Pt -k* p{fiq) p{ (q) d d q) 


1/2 


1 !/2 


< II Pt ll^ 2 ll fc * (Pt ~ Pt) lb < C^\\k * (jft - Pt ) 11 2 


n/ 


1 / 2 , 




(43) 


l+i 


From Lemma H.31 we know that as long as D{t ) < r N 2 , i.e. D*(t) < 1, 
the microscopic charge density is bounded as 


llrflloo <|B' , (2)l|lp[/,' , ]|U + B 2 (t) 


„“(V+ 2 ) n 2 / 


/•Ah 


<|H d (2)| sup ||p[// 

TVS N 

<|H d (2)|Co + l=:C p , 


+ 1 


(44) 


Note that this bound holds independent of N. Hence, we can use 
Loeper’s stability result, Proposition 14.11 in (1431) and get: 


\\k* (Pt ~Pt)h < [ max {ll^lloo, ||Hi 


flloo }YW 2 {plp{)<C2D{t). 


(45) 


Putting everything together and setting C\ := 2 C p Cl, we have 
d 


df 


D\t) < 2Ci(l V Vllog(rjv)l) D 2 (t ) 


or, after dividing by 2D[t) and multiplying both sides by r N 
d 


-(i+|) 


d t 


D*(t) < <7i(l v V|log(nv)|)£*(i). 


1/2 


19 

















By an application of Gronwall’s Lemma, we conclude that: 


D*(t ) < D*( 0) e l c ' 1 (v / l lo s( r Jv)l+ 1 )_ 

Finally, taking on the right hand side the inhmum over all 7To € 
/o), D*(0) becomes W*{p^ [Z\, /q) and we get for all t € T: 

W*(p?,f t N ) < W*(p^,f 0 )e tC ^\^(r N )\+i)_ (46) 

TTT / f\f X \ 

If there exists an e > 0 such that lim — Jf/dJ = 0, the right hand 

N^-oo rJ 

1 “I - — 

side converges to 0, so that, in particular, lim r N 2 W 2 {p^, f^) = 0. 

N—> oo 

□ 


To show convergence to solutions of the (unregularized) Vlasov-Poisson 
equation, we also require the following: 

Proposition 4.6. Let /o satisfy the assumptions of Proposition I SPA 
Let f t N and ft be the solution of the regularized, respectively the proper 
Vlasov-Poisson equation with initial data /o- Then: 

W 2 (f t N , ft) < r N e* Cl (Vl lo s(^)l+i). (47) 


Proof. Let p^ := p[ff N ] and pj := p[ft] be the charge density in¬ 
duced by f'l' 1 and ft, respectively. Let = (Q^,P t N ) the charac¬ 
teristic flow of f^ and tp t = ( Qt,Pt ) the characteristic flow of ft- We 
consider Tro(x,y) := fo(x)5(x — y) € II(/o,/o), which is already the 
optimal coupling yielding W.f (f t N , f t )\ t=0 = (Jo, ./o) = 0 and set 
7 r t = (ipt , ipt)#Tr 0 € LI (fl*, ft)- As above, we define 


D(t) : = 


J ((1 V \/|log(r/v)|) l^ 1 - y l \ + |x 2 - y 2 

R 6 xR 6 


2 I 1/2 

dv r t (x,y) 


and compute 


(48) 


^L» 2 (t) <2 J ((1 V Vllog^jv)!) | Q N (t, x) - Q(t, y)\ + | P N (t, x) - P(t, y)|) 

((1 V V|k>g(riv)|) \P N (t,x) - P(t, y)\ + \k*pf(Q N (x)) - k * p{ (QJy ))|) dn 0 (x,y) 

The proof proceeds analogous to Proposition 13.21 simplified by the 
fact that the charge densities remain bounded by assumption. The 
only noteworthy difference is in eq. m- Observing that k * p = 
k * p = k * (x N * X N * p), we use Lemma 14.41 to conclude: 


p t ) < W 2 {pp t ) + 2r N < W 2 (f t N , ft) + 2 r N < D(t) + 2 r N 

(49) 
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In total, we find: 


^ t D 2 {t) < 2 C 0 C L (1 V \/|log7 'jv|) D 2 {t ) + 2 C 0 D(t)(D(t) + 2 r N ) 
or 

< C\ (•y / |log r N | + 1 ) D(t) + 2 C 0 rAT 

with Ci > 2CqCl as defined in the previous proof. Using Gronwall’s 
inequality and the fact that D{ 0) = 0, we have 

W 2 (f t N J t ) < D(t) < r 7V e tCl (Vl 1 °g^l+ 1 ), 


from which the desired statement follows. 


□ 


4.3 Typicality 

To complete the proof of Theorem 13.31 it remains to show that the 
assumptions of Proposition 13.21 are satisfied for typical initial condi¬ 
tions, i.e. with probability approaching one as N tends to infinity. It 

is a classical result that if Z\, ..., Zjy are i.i.d. with law /, their empir- 

N 

ical density p, N [Z] = jj ^2 goes to / in probability. Establishing 

2=1 

quantitative bounds on large deviations (concentration estimates) is, 
however, a longstanding problem in probability theory with a vast 
amount of literature. To our knowledge, one of the first paper to ad¬ 
dress this question in the context of Wasserstein metrics was Bolley, 
Guillin, Villani, 2007 [3J. Subsequently, other authors have derived 
stronger concentration estimates, see, in particular, [2] and [5|. Very 
recently, great progress has been made in the paper of Fournier and 
Guillin, 2014, which considerably improves upon previous results, both 
in strength and generality [7]. We cite now their concentration esti¬ 
mates and apply them to conclude the proof of our main theorem. 

Theorem 4.7 (Fournier and Guillin). 

Let f be a probability measure on M n such that 3k > 2 p: 

M k (/) := J\x\ k df{x) < +oo. 

R fc 

Let {Zi)i = i n be a sample of independent variables, distributed ac- 

N 

cording to the law f and consider fi N [Z] := 22 ■ Then, for any 

2=1 
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e > 0 there exist constants c,C depending only on k,Mj e (f) and e such 
that for all N >1 and f > 0: 


with 


P 0 (w${p, N ,f) >f)< CN(Nf)- 1 ^ +Cts< 1 a(N,0 


a(N, £) 


'exp(— cNf 2 ) if P > n/2 

< exp(-cA^( In( 4i/g) ) 2 ) ifp = n /2 

^exp(— cN£ k / p ) if p £ [l,n/2). 


(50) 


We now apply this result to conclude the proof of our main theorem 
and establish an upper bound on rjv. 


Proof of Theorem 1 3.3\ Let rjy > N s and e > 0. Let A C be 
the (IV-dependent) set dehned by 

Z€A ^ W 2 (^[Z\, / 0 ) > r^ +£ . (51) 

We apply the previous in n = 2d dimensions with £ = _/\r-<H 2 + d + 2e ) < 
r' A ) 2 and the hnite moment assumption (1251) . We hnd: 


F 0 (A) < c(exp(-cNN~^ 2+d+2e)d ) + N 1 - 
Where the probability is dehned with respect to (& N fo- Choosing 


<5 


1 - e 

(2 + d + 2 e)d 


(52) 


we have 


P 0 (M) < C (exp(—clV 6 ) + N 1 - 1+&) -+ 0, N -+ oo. 


For the typical initial conditions Z € „4. c , it holds according to Propo¬ 
sition [3T2] that for all t < T. 

W*(p?,f t N ) < W*(p»,fo)e tC ^°s(r N )\ + i) 

< (IV v'|log(^)|)^ (1+ ^W 2 (^Jo)e tCl (Vl log ^)l +1 ) 

< (1 V V|log(r JV )|)r^e TC, 1 (v'l 1 °s(»-iv)|+i). ( 53 ) 


Observing that eV / l logrjv l = (e logr Jv) \/l lo snvl = (rjv) ^ /|logT ' ivl ) there 
exists ]Vo G N such that fl53l) < 1 for all N > Nq. More precisely, it 


suffices to choose Nq large enough that rjsr 0 < e 1 
hnd: 


2CiT+l .2 
e ' 


Then we 


l+l 


W*(p?jn < 1 =+ W 2 W,ft) < rJ'W'Qjfjn < r 


i+i 

'n 


(54) 
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Now we recall from Proposition 14.61 

W 2 (f t N J t ) < r N e «?i(Vli°e(™)l+i) 


which is smaller than \r l N e for N > iVo- We conclude the proof by 
noting that 


W 2 (^[Z\,f t )<W 2 (^[Z\ 


ft 


N ' + W 2 (f t N ,f t )<r 


1+i + -rlr £ < 


,1-e 


JV 


' AT 


for all Z E A c , TV > iVo and f E [0, T]. 


□ 


23 



References 


[1] N. Boers and P. Pickl. On mean field limits for dynamical systems. 
Journal of Statistical Physics, pages 1-16, 2015. 

[2] E. Boissard. Simple bounds for the convergence of empirical and 
occupation measures in 1-Wasserstein distance. Electronic Jour¬ 
nal of Probability, 16:2296-2333, 2011. 

[3] F. Bolley, A. Guillin, and C. Villani. Quantitative concentration 
inequalities for empirical measures on non-compact spaces. Prob¬ 
ability Theory and Related Fields, 137(3):541-593, 2007. 

[4] W. Braun and K. Hepp. The Vlasov dynamics and its fluctuations 
in the 1/N limit of interacting classical particles. Conun. Math. 
Phys., 56(2):101—113, 1977. 

[5] S. Dereich, M. Scheutzow, and R. Schottstedt. Constructive quan¬ 
tization: approximation by empirical measures. In Annales de 
I’Institut Henri Poincare, Probabilities et Statistiques, volume 49, 
pages 1183-1203. Institut Henri Poincare, 2013. 

[6] R. L. Dobrushin. Vlasov equations. Functional Analysis and Its 
Applications, 13(2) :115 123, 1979. 

[7] N. Fournier and A. Guillin. On the rate of convergence in wasser- 
stein distance of the empirical measure. Probability Theory and 
Related Fields, 162:1-32, 2014. 

[8] F. Golse. The mean-held limit for a regularized Vlasov-Maxwell 
dynamics. Commun. Math. Phys., 310(3):789-816, 2012. 

[9] F.A. Griinbaum. Propagation of chaos for the Boltzmann equa¬ 
tion. Arch. Rational Mech. Anal., 42:323-345, 1971. 

[10] M. Hauray and P.-E. Jabin. Particles approximations of Vlasov 
equations with singular forces : Propagation of chaos. To appear 
in Ann. Sci. Ec. Norm. Super, 2013. 

[11] E. Horst. Global strong solutions of Vlasov’s equation - neces¬ 
sary and sufficient conditions for their existence. Banach Center 
Publications, 19(1) :143—153, 1987. 

[12] E. Horst. On the asymptotic growth of the solutions of the Vlasov- 
Poisson system. Math. Methods Appl. Sci., 16:75-85, 1993. 

[13] M. Kac. Foundations of kinetic theory. In Proceedings of the Third 
Berkeley Symposium on Mathematical Statistics and Probability, 
1954-1955, volume vol III, pages 171-197. University of California 
Press, 1956. 

[14] D. Lazarovici and P. Pickl. A mean-held limit for the Vlasov- 
Poisson system. Preprint: arXiv:1502.04608, 2015. 


24 



[15] P.-L. Lions and B. Perthame. Propagation of moments and reg¬ 
ularity for the 3-dimensional Vlasov-Poisson system. Inventiones 
Mathematicae, 105:415-430, 1991. 

[16] G. Loeper. Uniqueness of the solution to the Vlasov-Poisson sys¬ 
tem with bounded density. J. Math. Pures Appl.. 86:68-79, 2006. 

[17] S. Mischler and M. Hauray. On Kac’s chaos and related problems. 
Journal of Functional Analysis , 266(10) :6055-6157, 2014. 

[18] H. Neunzert. An introduction to the nonlinear Boltzmann-Vlasov 
equation. In Carlo Cercignani, editor, Kinetic Theories and the 
Boltzmann Equation , volume 1048 of Lecture Notes in Mathemat¬ 
ics , pages 60-110. Springer Berlin Heidelberg, 1984. 

[19] H. Neunzert and J. Wick. Die Approximation der Losung 
von Integro-Differentialgleichungen durch endliche Punktmengen. 
In R. Ansorge and W. Tornig, editors, Numerische Behand- 
lung nichtlinearer Integrodifferential - und Differentialgleichun- 
gen, volume 395 of Lecture Notes in Mathematics , pages 275-290. 
Springer, Berlin, Heidelberg, 1974. 

[20] K. Pfaffelmoser. Global classical solutions of the Vlasov-Poisson 
system in three dimensions for general initial data. J. Differential 
Equations , 95(2):281-303, 1990. 

[21] G. Rein. Global weak solutions of the relativistic Vlasov-Maxwell 
system revisited. Comm, in Math. Sci., 2:145-158, 2004. 

[22] J. Schaeffer. Global existence of smooth solutions to the Vlasov- 
Poisson system in three dimensions. Comm. Partial Differential 
Equations, 16(8-9): 1313-1335, 1991. 

[23] H. Spohn. Large scale dynamics of interacting particles. Springer 
Berlin Heidelberg, 1991. 

[24] H. Spohn. Dynamics of Charged Particles and their Radiation 
Field. Cambridge University Press, 2004. 

[25] A.-S. Sznitman. Topics in propagation of chaos. In Ecole d’Ete de 
Probabilites de Saint-Flour XIX - 1989, volume 1464 of Lexture 
Notes in Mathematics, pages 165-251. Springer, Berlin, 1991. 

[26] S. Ukai and T. Okabe. On classical solutions in the large in time 
of two-dimensional Vlasov’s equation. Osaka J. Math., 15(2):245- 
261, 1978. 

[27] C. Villani. Optimal Transport Old and New, volume 338 of 
Grundlehren der mathematischen Wissenschaften. Springer, 
Berlin, 2009. 


25 



[28] S. Wollman. Global-in-time solutions of the two-dimensional 
Vlasov-Poisson systems. Communications on Pure and Applied 
Mathematics, 33(2):1T3 197, 1980. 


26 



